The influence of model building schemes and molecular dynamics sampling on QM-cluster models: the chorismate mutase case study.

Most QM-cluster models of enzymes are constructed based on X-ray crystal structures, which limits comparison to in vivo structure and mechanism. The active site of chorismate mutase from Bacillus subtilis and the enzymatic transformation of chorismate to prephenate is used as a case study to guide construction of QM-cluster models built first from the X-ray crystal structure, then from molecular dynamics (MD) simulation snapshots. The Residue Interaction Network ResidUe Selector (RINRUS) software toolkit, developed by our group to simplify and automate the construction of QM-cluster models, is expanded to handle MD to QM-cluster model workflows. Several options, some employing novel topological clustering from residue interaction network (RIN) information, are evaluated for generating conformational clustering from MD simulation. RINRUS then generates a statistical thermodynamic framework for QM-cluster modeling of the chorismate mutase mechanism via refining 250 MD frames with density functional theory (DFT). The 250 QM-cluster models sampled provide a mean ΔG‡ of 10.3 ± 2.6 kcal mol-1 compared to the experimental value of 15.4 kcal mol-1 at 25 °C. While the difference between theory and experiment is consequential, the level of theory used is modest and therefore "chemical" accuracy is unexpected. More important are the comparisons made between QM-cluster models designed from the X-ray crystal structure versus those from MD frames. The large variations in kinetic and thermodynamic properties arise from geometric changes in the ensemble of QM-cluster models, rather from the composition of the QM-cluster models or from the active site-solvent interface. The findings open the way for further quantitative and reproducible calibration in the field of computational enzymology using the model construction framework afforded with the RINRUS software toolkit.

Hydrophobic Clusters Regulate Surface Hydration Dynamics of Bacillus subtilis Lipase A.

The surface hydration diffusivity of Bacillus subtilis Lipase A (BSLA) has been characterized by low-field Overhauser dynamic nuclear polarization (ODNP) relaxometry using a series of spin-labeled constructs. Sites for spin-label incorporation were previously designed via an atomistic computational approach that screened for surface exposure, reflective of the surface hydration comparable to other proteins studied by this method, as well as minimal impact on protein function, dynamics, and structure of BSLA by excluding any surface site that participated in greater than 30% occupancy of a hydrogen bonding network within BSLA. Experimental ODNP relaxometry coupling factor results verify the overall surface hydration behavior for these BSLA spin-labeled sites similar to other globular proteins. Here, by plotting the ODNP parameters of relative diffusive water versus the relative bound water, we introduce an effective "phase-space" analysis, which provides a facile visual comparison of the ODNP parameters of various biomolecular systems studied to date. We find notable differences when comparing BSLA to other systems, as well as when comparing different clusters on the surface of BSLA. Specifically, we find a grouping of sites that correspond to the spin-label surface location within the two main hydrophobic core clusters of the branched aliphatic amino acids isoleucine, leucine, and valine cores observed in the BSLA crystal structure. The results imply that hydrophobic clustering may dictate local surface hydration properties, perhaps through modulation of protein conformations and samplings of the unfolded states, providing insights into how the dynamics of the hydration shell is coupled to protein motion and fluctuations.

Progress in Identification of UDP-Glycosyltransferases for Ginsenoside Biosynthesis.

Ginsenosides, the primary pharmacologically active constituents of the Panax genus, have demonstrated a variety of medicinal properties, including anticardiovascular disease, cytotoxic, antiaging, and antidiabetes effects. However, the low concentration of ginsenosides in plants and the challenges associated with their extraction impede the advancement and application of ginsenosides. Heterologous biosynthesis represents a promising strategy for the targeted production of these natural active compounds. As representative triterpenoids, the biosynthetic pathway of the aglycone skeletons of ginsenosides has been successfully decoded. While the sugar moiety is vital for the structural diversity and pharmacological activity of ginsenosides, the mining of uridine diphosphate-dependent glycosyltransferases (UGTs) involved in ginsenoside biosynthesis has attracted a lot of attention and made great progress in recent years. In this paper, we summarize the identification and functional study of UGTs responsible for ginsenoside synthesis in both plants, such as Panax ginseng and Gynostemma pentaphyllum, and microorganisms including Bacillus subtilis and Saccharomyces cerevisiae. The UGT-related microbial cell factories for large-scale ginsenoside production are also mentioned. Additionally, we delve into strategies for UGT mining, particularly potential rapid screening or identification methods, providing insights and prospects. This review provides insights into the study of other unknown glycosyltransferases as candidate genetic elements for the heterologous biosynthesis of rare ginsenosides.

Characterization of a transitionally occupied state and thermal unfolding of domain 1.1 of σ A factor of RNA polymerase from Bacillus subtilis.

σ factors are essential parts of bacterial RNA polymerase (RNAP) as they allow to recognize promotor sequences and initiate transcription. Domain 1.1 of vegetative σ factors occupies the primary channel of RNAP and also prevents binding of the σ factor to promoter DNA alone. Here, we show that domain 1.1 of Bacillus subtilis σ A exists in more structurally distinct variants in dynamic equilibrium. The major conformation at room temperature is represented by a previously reported well-folded structure solved by nuclear magnetic resonance (NMR), but 4% of the protein molecules are present in a less thermodynamically favorable state. We show that this population increases with temperature and we predict its significant elevation at higher but still biologically relevant temperatures. We characterized the minor state of the domain 1.1 using specialized methods of NMR. We found that, in contrast to the major state, the detected minor state is partially unfolded. Its propensity to form secondary structure elements is especially decreased for the first and third α helices, while the second α helix and ß strand close to the C-terminus are more stable. We also analyzed thermal unfolding of the domain 1.1 and performed functional experiments with full length σ A and its shortened version lacking domain 1.1 ( σ A _ Δ 1.1 ). The results revealed that while full length σ A increases transcription activity of RNAP with increasing temperature, transcription with σ A _ Δ 1.1 remains constant. In summary, this study reveals conformational dynamics of domain 1.1 and provides a basis for studies of its interaction with RNAP and effects on transcription regulation.

Comparative transcriptome analysis of wild type and a pectate lyase mutant strain of Bacillus subtilis.

Transcriptome analysis for the original Bacillus subtilis K1 strain and UV mutagenic strain UW07 with high yield of pectate lyase was implemented with RNA-seq. The function of genes was annotated and metabolic pathways were classified to look for different expression genes and classify these genes into related metabolic pathways to reveal the high-yield mechanism of pectate lyase in UW07. The results showed that 397 genes were up-regulated and 617 genes were down-regulated compared with the original strain. The up-regulated genes were mainly involved in ABC transporters, two-component system, biosynthesis of amino acids, and carbon metabolism.

Intramolecular electron transfer from biopterin to FeII-O2 complex in nitric oxide synthases occurs at very different rates between bacterial and mammalian enzymes: Direct observation of a catalytically active intermediate.

Nitric oxide synthase (NOS) is a cytochrome P450-type mono­oxygenase that catalyzes the oxidation of L-arginine to nitric oxide. We previously observed that intramolecular electron transfer from biopterin to Fe2+-O2 in Deinococcus radiodurans NOS (DrNOS) using pulse radiolysis. However, the rate of electron transfer in DrNOS (2.2 × 103 s-1) contrasts with a reported corresponding rate (11 s-1) in a mammalian NOS determined using rapid freeze-quench (RFQ) EPR. We applied pulse radiolysis to Bacillus subtilis NOS (bsNOS) and to rat neural NOS oxygenase domain NOS (mNOS). Concurrently, RFQ EPR was used to trap a pterin radical during single-turnover enzyme reactions of the enzymes. By using the pulse radiolysis method, hydrated electrons (eaq-) reduced the heme iron of NOS enzymes. Subsequently, ferrous heme reacted with O2 to form a Fe2+-O2 intermediate. In the presence of pterin, the intermediate of bsNOS was found to convert to other intermediate in the time range of milliseconds. A similar process was determined to have occurred after pulse radiolysis of the pterin-bound mNOS, though the rate was much slower. The intermediates of all of the NOS enzymes further converted to the original ferric form in the time range of seconds. When using the RFQ method, pterin radicals were formed very rapidly in both DrNOS and bsNOS in the time range of milliseconds. In contrast, the pterin radical in mNOS was observed to form slowly, at a rate of ∼20 s-1.

Structural Insights into the Dimeric Form of Bacillus subtilis RNase Y Using NMR and AlphaFold.

RNase Y is a crucial component of genetic translation, acting as the key enzyme initiating mRNA decay in many Gram-positive bacteria. The N-terminal domain of Bacillus subtilis RNase Y (Nter-BsRNaseY) is thought to interact with various protein partners within a degradosome complex. Bioinformatics and biophysical analysis have previously shown that Nter-BsRNaseY, which is in equilibrium between a monomeric and a dimeric form, displays an elongated fold with a high content of α-helices. Using multidimensional heteronuclear NMR and AlphaFold models, here, we show that the Nter-BsRNaseY dimer is constituted of a long N-terminal parallel coiled-coil structure, linked by a turn to a C-terminal region composed of helices that display either a straight or bent conformation. The structural organization of the N-terminal domain is maintained within the AlphaFold model of the full-length RNase Y, with the turn allowing flexibility between the N- and C-terminal domains. The catalytic domain is globular, with two helices linking the KH and HD modules, followed by the C-terminal region. This latter region, with no function assigned up to now, is most likely involved in the dimerization of B. subtilis RNase Y together with the N-terminal coiled-coil structure.

Biodiversity of N-acyl homoserine lactonase (aiiA) gene from Bacillus subtilis.

Microorganisms rely on the benefit of using chemical signals called autoinducers (AIs) as a connection matter in term of population, this mechanism is known as quorum sensing (QS). Quorum sensing is responsible for formation of biofilm together with virulence in bacteria. The most known QS molecule is N-acyl homoserine lactones (AHLs). A lot of degrading enzymes including lactonases that open the AHL ring and acylases that breakdown its acyl side chain can degrade or inactivate AHL. Due to similarity in lactone ring structure among AHLs it is susceptible to most of lactonases. Bacillus species are among the most promising bacteria producing AHL-lactonase. The aim of the work is to identify and study the diversity of the AHL-Lactonase gene among different Bacillus subtilis as a promising Quorum Quenching (QQ) strategy to prevent bacterial infections and biofilm formation. The AHL-lactonase (aiiA) gene of 64 B. subtilis isolates was amplified and sequenced followed by multiple sequence alignment of the translated amino acid sequences, homology modeling and docking study. An expected PCR product of about 750 base pair was detected in 22 B. subtilis isolates, and the results revealed that the isolates' sequences showed identity ranged between 97.61% and 99.47% with those in the NCBI GenBank database with 100% query coverage and 0.0 E-value. In addition, the results revealed high level of identity between many aiiA gene sequences of our isolates as they were closely related to the same sequences to many sequences of the NCBI GenBank database. The alignment of the amino acid sequences from the 22 B. subtilis isolates indicated that 84.4% of the amino acid residues were conserved between the aligned sequences. Docking of the co-crystalized ligand to wildtype and H109Y mutated protein showed a significant reduction of docking score for the mutated protein. This result indicate that this mutation might affect recognition or at least kinetics of these enzymes and hence their roles in quorum-quenching.

Protein-protein interactions enhance the thermal resilience of SpyRing-cyclized enzymes: A molecular dynamic simulation study.

Recently a technique based on the interaction between adhesion proteins extracted from Streptococcus pyogenes, known as SpyRing, has been widely used to improve the thermal resilience of enzymes, the assembly of biostructures, cancer cell recognition and other fields. It was believed that the covalent cyclization of protein skeleton caused by SpyRing reduces the conformational entropy of biological structure and improves its rigidity, thus improving the thermal resilience of the target enzyme. However, the effects of SpyTag/ SpyCatcher interaction with this enzyme are poorly understood, and their regulation of enzyme properties remains unclear. Here, for simplicity, we took the single domain enzyme lichenase from Bacillus subtilis 168 as an example, studied the interface interactions in the SpyRing by molecular dynamics simulations, and examined the effects of the changes of electrostatic interaction and van der Waals interaction on the thermal resilience of target enzyme. The simulations showed that the interface between SpyTag/SpyCatcher and the target enzyme is different from that found by geometric matching method and highlighted key mutations at the interface that might have effect on the thermal resilience of the enzyme. Our calculations highlighted interfacial interactions between enzyme and SpyTag/SpyCatcher, which might be useful in rational designs of the SpyRing.

A hydrophilic microenvironment in the substrate-translocating groove of the YidC membrane insertase is essential for enzyme function.

The YidC family of proteins are membrane insertases that catalyze the translocation of the periplasmic domain of membrane proteins via a hydrophilic groove located within the inner leaflet of the membrane. All homologs have a strictly conserved, positively charged residue in the center of this groove. In Bacillus subtilis, the positively charged residue has been proposed to be essential for interacting with negatively charged residues of the substrate, supporting a hypothesis that YidC catalyzes insertion via an early-step electrostatic attraction mechanism. Here, we provide data suggesting that the positively charged residue is important not for its charge but for increasing the hydrophilicity of the groove. We found that the positively charged residue is dispensable for Escherichia coli YidC function when an adjacent residue at position 517 was hydrophilic or aromatic, but was essential when the adjacent residue was apolar. Additionally, solvent accessibility studies support the idea that the conserved positively charged residue functions to keep the top and middle of the groove sufficiently hydrated. Moreover, we demonstrate that both the E. coli and Streptococcus mutans YidC homologs are functional when the strictly conserved arginine is replaced with a negatively charged residue, provided proper stabilization from neighboring residues. These combined results show that the positively charged residue functions to maintain a hydrophilic microenvironment in the groove necessary for the insertase activity, rather than to form electrostatic interactions with the substrates.

Structural and Biochemical Analysis of the Furan Aldehyde Reductase YugJ from Bacillus subtilis.

NAD(H)/NADP(H)-dependent aldehyde/alcohol oxidoreductase (AAOR) participates in a wide range of physiologically important cellular processes by reducing aldehydes or oxidizing alcohols. Among AAOR substrates, furan aldehyde is highly toxic to microorganisms. To counteract the toxic effect of furan aldehyde, some bacteria have evolved AAOR that converts furan aldehyde into a less toxic alcohol. Based on biochemical and structural analyses, we identified Bacillus subtilis YugJ as an atypical AAOR that reduces furan aldehyde. YugJ displayed high substrate specificity toward 5-hydroxymethylfurfural (HMF), a furan aldehyde, in an NADPH- and Ni2+-dependent manner. YugJ folds into a two-domain structure consisting of a Rossmann-like domain and an α-helical domain. YugJ interacts with NADP and Ni2+ using the interdomain cleft of YugJ. A comparative analysis of three YugJ structures indicated that NADP(H) binding plays a key role in modulating the interdomain dynamics of YugJ. Noticeably, a nitrate ion was found in proximity to the nicotinamide ring of NADP in the YugJ structure, and the HMF-reducing activity of YugJ was inhibited by nitrate, providing insights into the substrate-binding mode of YugJ. These findings contribute to the characterization of the YugJ-mediated furan aldehyde reduction mechanism and to the rational design of improved furan aldehyde reductases for the biofuel industry.

The WalR-WalK Signaling Pathway Modulates the Activities of both CwlO and LytE through Control of the Peptidoglycan Deacetylase PdaC in Bacillus subtilis.

The WalR-WalK two component signaling system in Bacillus subtilis functions in the homeostatic control of the peptidoglycan (PG) hydrolases LytE and CwlO that are required for cell growth. When the activities of these enzymes are low, WalR activates transcription of lytE and cwlO and represses transcription of iseA, a secreted inhibitor of LytE. Conversely, when PG hydrolase activity is too high, WalR-dependent expression of lytE and cwlO is reduced and iseA is derepressed. In a screen for additional factors that regulate this signaling pathway, we discovered that overexpression of the membrane-anchored PG deacetylase PdaC increases WalR-dependent gene expression. We show that increased expression of PdaC, but not catalytic mutants, prevents cell wall cleavage by both LytE and CwlO, explaining the WalR activation. Importantly, the pdaC gene, like iseA, is repressed by active WalR. We propose that derepression of pdaC when PG hydrolase activity is too high results in modification of the membrane-proximal layers of the PG, protecting the wall from excessive cleavage by the membrane-tethered CwlO. Thus, the WalR-WalK system homeostatically controls the levels and activities of both elongation-specific cell wall hydrolases. IMPORTANCE Bacterial growth and division requires a delicate balance between the synthesis and remodeling of the cell wall exoskeleton. How bacteria regulate the potentially autolytic enzymes that remodel the cell wall peptidoglycan remains incompletely understood. Here, we provide evidence that the broadly conserved WalR-WalK two-component signaling system homeostatically controls both the levels and activities of two cell wall hydrolases that are critical for cell growth.

Creation of an Engineered Amide Synthetase Biocatalyst by the Rational Separation of a Two-Step Nitrile Synthetase.

The synthesis of amides through acid and amine coupling is one of the most commonly used reactions in medicinal chemistry, yet still requires atom-inefficient coupling reagents. There is a current demand to develop greener, biocatalytic approaches to amide bond formation. The nitrile synthetase (NS) enzymes are a small family of ATP-dependent enzymes which catalyse the transformation of a carboxylic acid into the corresponding nitrile via an amide intermediate. The Bacillus subtilis QueC (BsQueC) is an NS involved in the synthesis of 7-cyano-7-deazaguanine (CDG) natural products. Through sequence homology and structural analysis of BsQueC we identified three highly conserved residues, which could potentially play important roles in NS substrate binding and catalysis. Rational engineering led to the creation of a NS K163A/R204A biocatalyst that converts the CDG acid into the primary amide, but does not proceed to the nitrile. This study suggests that NSs could be further developed for coupling agent-free, amide-forming biocatalysts.

Discovery and initial characterization of YloC, a novel endoribonuclease in Bacillus subtilis.

The Bacillus subtilis genome is predicted to encode numerous ribonucleases, including four 3' exoribonucleases that have been characterized to some extent. A strain containing gene knockouts of all four known 3' exoribonucleases is viable, suggesting that one or more additional RNases remain to be discovered. A protein extract from the quadruple RNase mutant strain was fractionated and RNase activity was followed, resulting in the identification of an enzyme activity catalyzed by the YloC protein. YloC is an endoribonuclease and is a member of the highly conserved "YicC family" of proteins that is widespread in bacteria. YloC is a metal-dependent enzyme that catalyzes the cleavage of single-stranded RNA, preferentially at U residues, and exists in an oligomeric form, most likely a hexamer. As such, YloC shares some characteristics with the SARS-CoV Nsp15 endoribonuclease. While the in vivo function of YloC in B. subtilis is yet to be determined, YloC was found to act similarly to YicC in an Escherichia coli in vivo assay that assesses decay of the small RNA, RyhB. Thus, YloC may play a role in small RNA regulation.

Chemoproteomics profiling of surfactin-producing nonribosomal peptide synthetases in living bacterial cells.

Much of our current knowledge on nonribosomal peptide synthetases (NRPSs) is based on studies in which the full NRPS system or each protein domain is expressed in heterologous hosts. Consequently, methods to detect the endogenous activity of NRPSs, under natural cellular conditions, are needed for the study of NRPS cell biology. Here, we describe the in vivo activity-based protein profiling (ABPP) for endogenous NRPSs and its applications to the study of their activities in bacteria. Remarkably, in vitro and in vivo ABPP in the context of the surfactin producer Bacillus subtilis enabled the visualization, tracking, and imaging of an endogenous SrfAB-NRPS with remarkable selectivity and sensitivity. Furthermore, in vivo, ABPP allowed the discovery of the degradation processes of the endogenous SrfAB-NRPS in the context of its native producer bacteria. Overall, this study deepens our understanding of the properties of NRPSs that cannot be addressed by conventional methods.

Distinct Interaction Mechanism of RNA Polymerase and ResD at Proximal and Distal Subsites for Transcription Activation of Nitrite Reductase in Bacillus subtilis.

The ResD-ResE signal transduction system plays a pivotal role in anaerobic nitrate respiration in Bacillus subtilis. The nasD operon encoding nitrite reductase is essential for nitrate respiration and is tightly controlled by the ResD response regulator. To understand the mechanism of ResD-dependent transcription activation of the nasD operon, we explored ResD-RNA polymerase (RNAP), ResD-DNA, and RNAP-DNA interactions required for nasD transcription. Full transcriptional activation requires the upstream promoter region where five molecules of ResD bind. The distal ResD-binding subsite at -87 to -84 partially overlaps a sequence similar to the consensus distal subsite of the upstream (UP) element with which the Escherichia coli C-terminal domain of the α subunit (αCTD) of RNAP interacts to stimulate transcription. We propose that interaction between αCTD and ResD at the promoter-distal site is essential for stimulating nasD transcription. Although nasD has an extended -10 promoter, it lacks a reasonable -35 element. Genetic analysis and structural simulations predicted that the absence of the -35 element might be compensated by interactions between σA and αCTD and between αCTD and ResD at the promoter-proximal ResD-binding subsite. Thus, our work suggested that ResD participates in nasD transcription activation by binding to two αCTD subunits at the proximal and distal promoter sites, representing a unique configuration for transcription activation. IMPORTANCE A significant number of ResD-controlled genes have been identified, and transcription regulatory pathways in which ResD participates have emerged. Nevertheless, the mechanism of how ResD activates transcription of different genes in a nucleotide sequence-specific manner has been less explored. This study suggested that among the five ResD-binding subsites in the promoter of the nasD operon, the promoter-proximal and -distal ResD-binding subsites play important roles in nasD activation by adapting different modes of protein-protein and protein-DNA interactions. The finding of a new type of protein-promoter architecture provides insight into the understanding of transcription activation mechanisms controlled by transcription factors, including the ubiquitous response regulators of two-component regulatory systems, particularly in Gram-positive bacteria.

DisA Restrains the Processing and Cleavage of Reversed Replication Forks by the RuvAB-RecU Resolvasome.

DNA lesions that impede fork progression cause replisome stalling and threaten genome stability. Bacillus subtilis RecA, at a lesion-containing gap, interacts with and facilitates DisA pausing at these branched intermediates. Paused DisA suppresses its synthesis of the essential c-di-AMP messenger. The RuvAB-RecU resolvasome branch migrates and resolves formed Holliday junctions (HJ). We show that DisA prevents DNA degradation. DisA, which interacts with RuvB, binds branched structures, and reduces the RuvAB DNA-dependent ATPase activity. DisA pre-bound to HJ DNA limits RuvAB and RecU activities, but such inhibition does not occur if the RuvAB- or RecU-HJ DNA complexes are pre-formed. RuvAB or RecU pre-bound to HJ DNA strongly inhibits DisA-mediated synthesis of c-di-AMP, and indirectly blocks cell proliferation. We propose that DisA limits RuvAB-mediated fork remodeling and RecU-mediated HJ cleavage to provide time for damage removal and replication restart in order to preserve genome integrity.

Bacillus subtilis ANSB168 Producing d-alanyl-d-alanine Carboxypeptidase Could Alleviate the Immune Injury and Inflammation Induced by Ochratoxin A.

Ochratoxin A (OTA) is toxic to animals and threatens food safety through residues in animal tissues. A novel degrading strain Bacillus subtilis ANSB168 was isolated and further investigated. We cloned d-alanyl-d-alanine carboxypeptidase DacA and DacB from ANSB168 and over-expressed them in Escherichia coli Rosetta (DE3). Then, we characterized the OTA degradation mechanism of DacA and DacB, which was degrading OTA into OTα. A total of 45 laying hens were divided into three equal groups. The control group was fed basal feed, and other groups were administered with OTA (250 µg/kg of feed). A freeze-dried culture powder of ANSB168 (3 × 107 CFU/g, 2 kg/T of feed) was added to one of the OTA-fed groups for 28 days from day one of the experiment. We found that OTA significantly damaged the kidney and liver, inducing inflammation and activating the humoral immune system, causing oxidative stress in the layers. The ANSB168 bioproduct was able to alleviate OTA-induced kidney and liver damage, relieving OTA-induced inflammation and oxidative stress. Overall, DacA and DacB derived from ANSB168 degraded OTA into OTα, while the ANSB168 bioproduct was able to alleviate damages induced by OTA in laying hens.

Regulate the hydrophobic motif to enhance the non-classical secretory expression of Pullulanase PulA in Bacillus subtilis.

Bacillus subtilis has been widely used as a prokaryotic host for the secretory expression of heterologous proteins. In this study, a pullulanase (PulA) from Anoxybacillus sp. LM18-11 was firstly identified to be expressed in Bacillus subtilis 1A751 through non-classical secretion pathway. Results showed that both the N- and C-terminal regions of PulA were essential for its soluble expression. To explore its specific structural basis of secretion in B. subtilis, we revealed a hydrophobic motif A501-H507 which is vital for the secretion of the whole protein of PulA. Through a series of site-specific mutagenesis, the triple-sites mutants R503E/I506E/H507E and R503E/I506Y/H507E showed the highest extracellular activity (160.07 U/mL) and total activity (243.37 U/mL) which was 1.71 times and 1.55 times higher than those of PulA. The highest secretion rate of mutant I506E/H507E was more than 50% which was 34.72% higher comparing with that of PulA. The glutamic acid substitution on these three key surface sites which decreased the surface hydrophobicity of that region was confirmed to be beneficial to improve the secretory expression of PulA. This novel discovery for the secretory expression of PulA in B. subtilis would make a new perspective on regulating a kind of non-classical secretion in B. subtilis.